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LOCAL THERMAL EQUILIBRIUM FOR CERTAIN STOCHASTIC 
MODELS OF HEAT TRANSPORT 

YAO LI, PETER NANDORI, AND LAI-SANG YOUNG 


Abstract. This paper is about nonequilibrium steady states (NESS) of a class 
of stochastic models in which particles exchange energy with their “local environ¬ 
ments” rather than directly with one another. The physical domain of the system 
can be a bounded region of for any d > 1. We assume that the temperature at 
the boundary of the domain is prescribed and is nonconstant, so that the system 
is forced out of equilibrium. Our main result is local thermal equilibrium in the 
infinite volume limit. In the Hamiltonian context, this would mean that at any lo¬ 
cation X in the domain, local marginal distributions of NESS tend to a probability 
with density permitting one to define the local temperature at x to be 

j3{x)~^. We prove also that in the infinite volume limit, the mean energy profile 
of NESS satisfies Laplace’s equation for the prescribed boundary condition. Our 
method of proof is duality: by reversing the sample paths of particle movements, 
we convert the problem of studying local marginal energy distributions at x to that 
of joint hitting distributions of certain random walks starting from x, and prove 
that the walks in question become increasingly independent as system size tends to 
infinity. 


1. Introduction 

This paper attempts to address, using highly idealized models, two of the major 
challenges in nonequilibrium statistical mechanics: One is the derivation of the Fourier 
Law, equivalently the heat equation, from microscopic principles. The other has to 
do with local thermal equilibrium (LTE) for systems that are driven out of, and 
possibly far from, equilibrium. Both of these topics are of fundamental importance, 
and have been the focus of much research activity in recent years (see below), yet no 
satisfactory general theory has been proposed. In this paper we study the invariant 
measures corresponding to nonequilibrium steady states (NESS) of a specific class 
of particle systems with stochastic interactions. The models we consider are simple 
enough to be amenable to rigorous analysis, yet not overly specialized, so they may 
offer insight into more general situations. A feature that distinguishes our paper from 
a number of previous works is that the physical space of our models can be for 
any d > 1. Model behavior depends on d, necessitating different arguments in the 
proofs, but our results are valid for all d > 1. 

We begin with a rough model description; see Sect. 2.1 for more detail. For any 
d > 1, let D C be a bounded domain with smooth boundary, and let T be 
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a prescribed temperature function on dT>. We consider ~ particles performing 
independent random walks on Vi = fl LT) where L 3> 1 is a real number and LT) 
is the dilation of V. These particles do not interact directly with one another, but 
only via their “local environments”, symbolized by a collection of random variables 
each representing the energy at a site in T>l. More precisely, each particle carries 
with it an energy. As it moves about, it exchanges energy with each of the sites 
it visits, and when it reaches the boundary of T>l, it abandons the energy it was 
carrying, replacing it by an energy drawn randomly from the “bath distribution” at 
the corresponding point in dV. 

These stochastic models are modihcations of the 1-D mechanical chains studied in 
[EY06j and their 2-D generalizations in |LY10j . In these mechanical models, energy 
transport occurs via particle-disk interactions, an idea borrowed from |LLM03j . More 
precisely, there is an array of rotating disks evenly spaced in the domain. Particles 
do not interact with each other directly; they exchange energy with these disks upon 
collision. Our site energies are an abstraction of the kinetic energies of these disks, or 
the “tank energies” in the stochastic models in |EY06] . Further simplihcations have 
been introduced in the present models to make the analysis feasible. 

In our models, when the prescribed temperature function is constant on dT>, i.e. 
T = Tq for some Tq G M"*", it is easy to see that the unique invariant probability distri¬ 
bution is a product measure: particle numbers are independent and Poissonian, and 
particle and site energies are independent and exponentially distributed with mean 
Tq. Let us refer to such a distribution as “the equilibrium distribution at temperature 
To”. For a nonconstant T (all that we require is that it be a continuous function), 
the system is forced out of equilibrium. It is not hard to see that there is still a 
unique NESS to which all initial states converge. Our two main results assert that 
the following hold in the inhnite volume limit: 

(1) mean energy prohles with respect to NESS when scaled back to P 
converge to the unique solution u{x) of 

Am = 0 on "D , u\qx> = T ; 

(2) given any a; G T, for sites located near xL, marginal distributions of 
the NESS tend to the equilibrium distribution at temperature u{x). 

These and other results are formulated precisely in Sect. 2.2. 

Our method of proof is duality, and the dual used here is similar to that in |RY07j , 
which in turn borrowed its main idea from [KMP82j . We differ from these earlier 
works in that the KMP-model is not a particle system, and our results are proved for 
all d > 1. As in [KMP82j . our “dual process” (duality is with respect to a function, 
to be precise) keeps track of movements of certain discrete objects we call “packets” 
in this paper. Reasoning naively, marginal energy distributions at a site v G Tl is 
determined by what particles bring to this site from the bath, so of interest are the 
points of origin of these energies. The idea is that to identify these points of origin, 
we can place some packets at site v, to be carried around by particles in a manner 
analogous to the way energy is transported, run the particle trajectories “backwards”, 
and look at the hitting distributions of the packets on the boundary of As it turns 
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out, packet movements in this process are effectively independent random walks when 
their trajectories do not meet. From this, one observes that the problem is simpler in 
d > 3, where independent random walkers tend not to meet, and that for d = 1, where 
they meet often, exchangeability, which was used in [KMP82j . is entirely natural. 

Related works. We hnish with a discussion of some results that are related to 
the various aspects of the present paper; needless to say, this is very far from a 
complete list. Some general references on interacting particle systems are [L851IS911 
KLlSl. A sample of the articles on NESS are |RLL671 IKLS841 IEPR99i iRTiM IDLSn2i 
ID07|. IFLM11|. IKY13|. ICE141ILY14] . These and other works treat the existence and 
uniqueness of NESS, correlation decay, fluctuations etc. in stationary states. Two 
informative reviews on the Fourier Law are IBLRBOOj and |LLPn3j . Mathematical 
papers on the Fourier Law proven under various assumptions include |BO05t IBLL04 
IBK071 IRlH IDL11|. IL012|. IDN14j . There are many papers in the physics literature; 
see the references in |LLP03j and e.g. |GG08j . 

Turning now to LTE, for the physics of this phenomenon, see |GM62j . Our liter¬ 
ature search turned up relatively few rigorous results on LTE: An early important 
example is |KMP82j . see also |RY07j and |OR13j : a systematic study of certain gra¬ 
dient models in [ELS901 IELS91] : and the survey article |CGGR13] . It has been 
noted that a number of systems with integrable dynamics do not have LTE, see e.g. 
[DD991 IBLYlOl IRY12] . In [EY06j and |LY10] , the authors derived formulas for macro¬ 
scopic observations such as energy and particle density profiles for both mechanical 
and stochastic models assuming LTE, and provided numerical validation for their 
derivations. 

As for duality, two classical references are |S70j and |L85j . Two early influential 
papers in the use of this technique are |KMP82j and |PS83j : later works include 
e.g. jGKR07|. IBCS12] . A recent survey of duality from the probabilistic viewpoint is 
[JK14j . and a nice survey on the use of duality in statistical mechanics problems is 
[CGGRl^ : this survey contains many other references. 

2. Model and Results 

2.1. Model description. We hx a dimension d > 1, and let P C be either 
a bounded open rectangle or a bounded, connected open set for which dV is a 
submanifold. A continuous function T : dV —)■ R+ to be thought of as temperature is 
prescribed. For L S> 1, the physical domain of our system is 

VL = LVnz‘^, 

and each lattice point in Vl is referred to as a site. The hath is located at 

Bl = &Dl\Vl 

where dT>i = {u G : u has a neighbor in Vi and a neighbor outside of T>i}. 
Throughout this article, | ■ | denotes the cardinality of a hnite set. For any v G M'^, 
(u) denotes its closest point in and ||u|| is the Euclidean norm. 
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We consider a Markov process Xt = xj;^^ with random variables 
^ (rji,V m), (Xi, ...,Xm)) e x x 

where denotes the energy at site v G M = M{L) G is the nnmber of 
particles in the system, and rji and Xi are the energy and location of particle i. The 
inhnitesimal generator of this process has the form 

{Gf){x) = {GJ){x) + {G 2 f){x), 

where Gi and G 2 describe respectively interactions within Vl and with the bath: 
Each particle carries an exponential clock which rings at rate 1 independently of 
the clocks of other particles. When its clock rings, the particle exchanges energy 
with the site at which it is located; immediately thereafter it jnmps to a neighboring 
site, choosing the 2d nearest neighbors with eqnal probability. The Gi part of the 
generator describes the action when the neighboring site chosen is in 

M 1 

{Gif){x) = Y. Yd 

k=l w&VL-.\Xk-w\=l'^^ 

f{{Ov&VL,Vu-,Vk-u -p){^x,+Vk),Vk+i,-,VM, 

Xi, ..., Xfc_i, w, Xfc+ 1 ,..., Xm) — fix) , 

where 

/1 \ y' _ / fv X ^ Xj, 

~ 1 Pifx^ + Vk) if = Xfc. 

If the particle jnmps to a site v G Bl, then its energy is npdated according to the 
temperatnre at v, and it is retnrned immediately to its original site. More precisely, 
we extend T to a neighborhood of &D, and dehne 

1 /* 1 /*oo 

(G./)(x) = Em E II 

k=l weBL.\Xk-'w\ = l''^ 

/((Ot;eni,hD--,hfc-Dh',hfc+D-,hM,Xi,...,XM) - /(x) , 

where is given by ([I]) and /3(f) = T(f )-h 
This completes the dehnition of onr model. 

Remark on distinguishable vs indistinguishable particles. As dehned, the particles 
in Xt are named and distingnishable. Since onr resnlts pertain to inhnite-volnme 
limits, it is natnral to work with models with indistinguishable particles. For each 
L, snch a model can be obtained from x[G via the following identihcation: for 
x,y E X X T>ff, we let x^ y if 

X = ii^y)y^VL,iVl,-,VM),iXi,...,XM)) 
and y ) (i/(t(i)) • ) i7t7(Ar))) (Xa-(i),..., XCT(Ar))) 
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where a is a permutation of the set {1, 2,..., M}. It is easy to check that the quo¬ 
tient process is well dehned and corresponds to with indistinguishable 

particles. As the desired results for are deduced easily from those for xj:^\ 

we will, for the most part, be working with xj:^\ 

2.2. Statement of results. We begin with a result on the existence and uniqueness 
of invariant measure in the equilibrium case, i.e., when the prescribed bath tempera¬ 
ture T is constant. 

Proposition 1. Let d G Z’*' and V G be as above, and let L be such that 
restricted to Vl is a connected graph. If the function T is constant, then with the 
notation /3 = l/T, 

M M / 

v&V]^ j=l k=l V ^ v&Vl 

is the unique invariant probability measure of the process xj:^'^. 

We are primarily interested in out-of-equilibrium settings dehned by non-constant 
bath temperatures. 

Proposition 2. Let d,II and T be as in Sect. 2.1. We assume L is such that T/' 
restricted to Vl is a connected graph. Then the process has a unique invariant 
probability measure . Furthermore, the distribution of xj:^'^ converges to as 
t —)■ oo for any initial distribution of Xq^^ . 

The proofs of Propositions [T] and [2] are straightforward; thus we mention only the 
ideas and leave details to the reader. For Propositions [H one can check by a direct 
computation that is invariant; uniqueness follows from Doeblin’s condition. Once 
Proposition [1] is established, tightness (at “cx)” and at “0”) for Proposition [2] can be 
proved as follows: Given xj:^'^ dehned by a continuous bath temperature function 
T\dv, we consider two equilibrium processes corresponding to boundary conditions 
jnmax _ and = mix^sx>T{x). Let /rifliax and /ig^in be the invariant 

probabilities of these two processes respectively. Coupling xf^^ to the process dehned 
by in the natural way with Xg^^ = /rifliax, we see that the distribution of xf^^ 
is stochastically dominated by Likewise, coupling to the process dehned by 

jnmm, ggg distribution of xj:^'^ is stochastically bounded from below by 

(L) 

/^e,min* 

In Theorems 1-4 below, it is assumed that d,T>,T and L are as in Sect. 2.1. For 
given L, we let denote expectation with respect to the nonequilibrium steady 

state Our hrst result is on the prohle of mean site energies. 

Theorem 1 (Mean site energy profiles). For each x &V, 

lim = u{x) 
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where u is the unique solution of the equation 

(2) Am = 0 on D , M|aD = T. 

Next we proceed to the definition of LTE for site energies. For x G 'D and a finite 
set S C we denote by the projection of to the coordinates (^{a:L>+D)'!;es 
and identify with a measure on with coordinates Given M{L) G Z+, 

we say the site energies of approach local thermodynamic equilibrium (LTE) 

as L —>■ CX3 if for every x G and every finite set S' C Z'^, 

^ as L cx) , 

seS 

where ^ stands for weak convergence and /3{x) = u{x)~^, u being the function in 
Theorem [TJ 

Theorem 2 (LTE for site energies). For any choice of M{L) G Z"*", site energies 
of approach LTE as L —)■ oo . 

In the case where the number of particles tends to a fixed positive density as 
L ^ oo, i.e. , 

(3) — avol{V) as L —)■ oo 

Jj^ 

for some constant a > 0, a more complete notion of LTE should include not only 
distributions of site energies but also those of particle energies. In preparation for the 
formal statement, we introduce the following notation: 

Let be the steady state distribution of the process x]^'^ /^ (see the Remark 

at the end of Sect. 2.1). Given x G P, S' = {mi,...,Ms} C Z'^, and non-negative 
integers iLi, ...,Ks, we consider the conditional probability 

(/i^-^V-) I particles at site (xL) + Vj = Kj , j = 1,.. ., s} 
and project this measure to the site and particle energy coordinates on (xL) + S. The 
resulting measure, be viewed as a measure on 

X (M+^)/....), where the relation in (M;^)/.... is defined by x ~ i/ if x = 
{xi, ...,xk), y = (xo-(i),..., To-(ii')), and cr is a permutation of the set We 

are concerned with the limit of Ki Ks as L —?• oo. 

Theorem 3 (LTE for systems with positive density of particles). Suppose 
M{L)/L'^ —)■ avol(P) for some a > 0. Then the system approaches LTE as L ^ oo 
in the sense of both site and particle energy distributions. That is to say, for every 
X & V, S = {mi,...,Ms} C Z*^, and non-negative integers Ki,...,Ks, we have the 
following limiting distributions as L ^ oo: 

(1) Let be the number of particles at site (xL) + vj, j = 1,2,... ,s, seen as 
a random variable with respect to Then ... ki^'^ tend in distribution 

to independent Foisson random variables with mean a. 
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(2) The measures ^ converge weakly to 

(4) {8, X {8^^)/J) as L ^ oo , 

where 8^ is the exponential distribution on M+ with parameter (3 = u{x)~^, 8^ 
is the product of K copies of 8^, and ~ is the usual identification in 

The notion of LTE considered so far describes marginal energy distributions in 
regions of microscopic sizes. These results can be extended to a version of LTE at 
mesoscopic scales: 

Theorem 4 (LTE at mesoscopic scales). For any t) E (0,1), the results in Theo¬ 
rems 2 and 3 remain valid if we replace 

“ {{xL) + where S G is any finite set ” 

by “ {{xL + L'^v)).ues where S GMf is any finite set ” . 


3. Preliminaries from probability theory 

We collect in this section some facts from probability theory that will be used. All 
the results cited are known, possibly with the exception of Proposition [3l 


3.1. Random walks. Here, we formulate some basic lemmas about random walks. 
We will use the terminology random walk for any Markov chain of the form 


Sn = 


k=0 


where the fkS are independent, identically distributed random variables with values 
in {d is the dimension of the random walk). The special case where is supported 
on the origin and its 2d nearest neighbors is called a nearest neighbor random walk, 
and the case where is uniformly distributed on the 2d nearest neighbors of the 
origin is called a simple symmetric random walk (SSRW). The following statement is 
arguably the most important property of random walks with finite variance. 


Lemma 1 (Invariance principle). Consider a d-dimensional random walk Sn = X]fc=o 
where has zero expectation and finite covariance matrix S, and let Wn = (lTn(f))tg[o,i] 
be the random process defined by ITn(^) = ^ and linear interpolations between ^ and 

Then as n ^ oo, Wn converges weakly to the d-dimensional Brownian motion 
on [0,1] with covariance matrix S. 

Lemma 2 (Harmonic measure). Let x E F GMf where F is as above, and let r be 
the first hitting time of dF for a Brownian motion Bf starting from x. Then 

E(r(R^)) = u{x) 


where u is given by 
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The next result is a combination of the last two, together with a small perturbation 
in starting location. A version of this result likely exists in the literature, but we are 
unable to locate a reference. Since our proofs rely heavily on Proposition [3l we have 
included its proof in the Appendix. 


Proposition 3. Let x ^ V G be as above, and let e > 0 be given. Then there exist 
h > 0 and Lq such that the following holds true for all L > Lq: Let Sn be a SSRW on 
1/ with ll^o — xL\\ < 5L, and let r be the smallest n such that Sn G Then 


(5) 



u{x) 


< e. 


We will also use the following estimate on moderate deviations. It is a consequence 
of e.g. Theorem 1 in |P75j . Section VIII.2: 

Lemma 3. Let fn be bounded i.i.d real-valued random variables with variance cx^ > 0. 
There is a constant Ci such that for any n > 1 and for any R with 1 < R < nf!'^, 

F(|X^^fc-nE(0| 

\ k=l 

The reason for the present review is that instead of studying directly, we 

will transform the problem into one involving certain stochastic processes in which a 
hnite number of walkers perform SSRW on These walks are independent when 
the walkers are at distinct sites, but when they meet, there is a tendency for them to 
stick together for some random time. We will need to show that in terms of hitting 
distributions such as those in Lemma [21 the situation in the L ^ oo limit is as though 
the walks were independent. This is clearly related to the question of how often two 
walkers meet, a property well known to be dimension dependent. 

(1) d = 1. The local time of a SSRW in dimension 1 at the origin up to time n is 
~ y/n. More precisely, the local time up to n, rescaled by y/n converges weakly 
to the absolute value of the standard normal distribution (see |CH49] ). Since 
the typical time needed to leave the interval [—L, L] is 0{Lf) by the invariance 
principle, two independent walkers will meet ~ L times before leaving the 
interval [—L,L\. 

(2) d = 2. The local time of a SSRW in dimension 2 at the origin up to time n 
is ~ log n, so that two independent walkers meet ~ log L times before leaving 
the domain Vl. 

(3) d > 2. SSRWs in dimensions d > 2 are transient, meaning two independent 
walkers only meet finitely many times. 

These observations have prompted us to proceed as follows: We will hrst treat the 
d = 2 case, by comparing the process in question to independent walkers. Once that 
is done, we will observe that a simplified version of the argument gives immediately 
results for d > 2. Dimension one is treated differently: The large number of encounters 
makes it difficult to compare the stochastic processes above to independent walkers. 


< cie ^ 
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Instead, we make use of the meetings of the walkers to show that their identities can 
be switched; see Sect. 3.3. 

We collect below two other dimension-dependent facts about random walks that 
will be used in the sequel. 

Lemma 4 (d = 2). Let "D C he above. We fix /3 G (0,1), C < oo and e > 0, and 
let Wn he a SSRW on . Then 

(6) Pfmin{||iy„|| >L} <min{W„ = 0} | HWoH G > fi - e 

for all L large enough. 

It is a well known fact in the probabilistic literature that the left hand side of (]6|) 
converges to /? as L —)■ oo (see for instance Proposition 1.6.7 in |L91] b This implies 
Lemma 01 

For completeness, we provide a heuristic justification for this cited result: Note 
that Wn converges to a Brownian motion after rescaling. Since the logarithm of a 
planar Brownian motion is a martingale, it is not difficult to deduce the statement 
that if Bt is a planar Brownian motion with Bq 0, then the probability that ||Pt|| 
reaches 2||Po|| before reaching ||Po||/2 is 1/2. Consequently, 

(7) P ( ||lTn|| reaches 2L^ before reaching L^/2) ~ 1/2. 

For k E Z, let 

C^k = {v eZ^ : ||n|| < 2'^} , 

and define the random variables tj , kj in the following way: 

ti = min{3/c : Wn E dC 2 k} 

n>0 

and ki is such that Wt^ E dC 2 ki. Inductively, we define 
tj = min {Wn E or Wn E 

and kj such that Wf. E dC^^k^. Similarly to ([7]), we see that kj ^ log 2 ||fFq|| is 
approximately a one dimensional SSRW for L large with starting position fi log 2 L. 
A simple computation (often referred to as gambler’s ruin) gives that 

P (min{log 2 ||1W„|| = \ 0 g 2 L} < min{log 2 \\WtJ\ = 0}) fi . 

The following is a well known property of high dimensional random walks. 

Lemma 5 (d > 2). In dimensions d > 2, any non-degenerate random walk is tran¬ 
sient, i.e., ||S'„|| —)■ 00 as n —)■ 00 with probability 1. 

3.2. Moments of exponential random variables. Let A G M and s G Z'^. The 

moments of s independent exponential random variables Xi,..., Xg with parameter 
A are given by 

m(n„ .... n,) = E(AT‘.. . X;-) = n ^ , 

i=l 


(8) 
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Conversely, if Xi,... are such that their joint moments are given by ([8]) for all 
(rii,..., Hg) G {0, 1 , 2 ,... }^, then they are independent exponential random variables 
with parameter A. See for instance [Z83j . 

3.3. Exchangeability. We also review very briefly the notion of exchangeable ran¬ 
dom variables, which will be used in the sequel. 

Definition 1. The infinite sequence of random variables Xi, X 2 , ... is called exchange¬ 
able if for any finite N and any permutation a 0 / {1, 2,..., X}, the random vectors 
(Xi,...,X^r) and (Xo-(i),..., Xo-(Ar)) have the same distributions. 

The following result is well known; it was proved by de Finetti in [D31] : 

Theorem 5 (de Finetti’s Theorem). If Xi, X 2 , ... is a sequence of {0,1} -valued 
exchangeable random variables, then there exists a distribution function F on [0,1] 
such that for all n and all Xt G {0,1}, 

P{X,=Xi,...,Xn=Xn)= [ 

Jo 

4. Duality 

In this section we introduce another process Yj and a function F, and show that 
Xt and Yt are dual with respect to the function F in a sense to be made precise. We 
explain also how to leverage duality to prove some of the asserted results in Sect. 2.2. 

4.1. Motivation and definition of a “dual process”. The idea is as follows: For 
the process Xt = xj^\ consider the marginal energy distributions of at n G T>l- 
These distributions reflect what the particles bring to site v from the bath, and that 
in turn is reflected in which parts of jBl a particle visits prior to its arrival at site 
V — though in reality things are a bit more complicated: particles interact with all 
the sites they pass through on their way to v. Nevertheless, accepting this simplihed 
picture for the moment, we can reverse the trajectories of the particles, and think of 
them as carrying certain “packets” from site v to the bath. The locations in at 
which these packets are deposited will then tell us from which parts oi Bl energies 
were drawn in the original process, thereby revealing the composition of the steady 
state distribution of Xt at v. An advantage of studying the reverse process is that 
it is reminiscent of hitting distributions of Brownian motion, or harmonic measures 
on dD. These are the ideas behind the duality formulation discussed in this section. 

For T> <Z and L as in Sect. 2.1, we now introduce a Markov process 
designed to carry packets from the sites in Vl to Bl. This process also involves 
M = M{L) particles. Let N denote the set of non-negative integers. The variables in 
are 

n= {{nv)v(^VL, {hy)v(zB^,ni,...,nM,Yi,...,YM) G x x x . 
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Here, is to be interpreted as the number of packets at site v G 'T’l, the number 
of packets at v E Bl, hj the number of packets carried by particle j, and Yj the 
position of particle j. (We distinguish between packets that have been dropped off 
at a site and packets that are carried by particles: Uy counts packets that have been 
dropped off at site v, and does not include packets carried by particle j even when 
Yj = V.) The generator of the process Yt is given by 

{Af){n) = {AJ){n) + {A 2 f){n), 


where Ai corresponds to movements inside Vl and A 2 is the part describing the 
interaction with the boundary. Formally, we have 


iAif){n) 


E 


1 


E 

z&VL-.\Yk-z\=l 


1 

hfc + + 1 


fik+riz 

E IT- ■ •) - /(Dl 


q=0 


where the quantity inside the parenthesis is 




^/c—1: (ll'/c T 2 : < 7 ) y R/c+1 y •••y'^My : • • •: ^k—\y ^y y •" y 


with 

(9) 


{ Uy li V ^ z 

y q it V = z. 


Note the difference between ([T]) and ([9]): In when a particle’s clock rings, it pools 
its energy together with that at the site at which it is located, and carries a random 
fraction of it as it jumps to a new site. In Yt, the order is reversed: when a particle’s 
clock rings, it jumps to a new site, then pools together the packets it is carrying with 
those at the new site, and takes a random fraction, assuming this new site 2 ; is inside 
Vl. The second half of the generator treats the case where z E 3^: 


(^2/)(n) 


E 


1 


E 


1 

life + 1 




q=0 


where the quantity inside the parenthesis is 

ij^v)'v&T>L y ; hi, ..., flk—l, R-fc q, hfc-|_i, ..., , Yi, ..., Y^ , 


with 




fly \i V ^ z 
fly + hk it V = z 


and 


Hy if T 7^ Yfc 

q if T = Yfc- 


That is to say, if particle k jumps from site v E T>l to site 2 ; E El, then the following 
occurs instantaneously: it drops off all of the packets it is carrying at site z, returns to 
site V, and takes a random fraction of the packets located at site v. Once a packet is 
dropped off in Bl, it will remain there permanently, so that as time tends to infinity 
there will be no packets left in 


This completes the definition of the process 


v-T) 
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4.2. Proof of pathwise duality. The function with respect to which duality will 
be proved is 


M 


v & Vl •’=1 ^ 


i=i 


v&Bl 


where x, and rjj are as in the dehnition of Xt and the rest are from the dehnition 
of Yt- Notice that F does not depend on the positions of the particles in either Xt 
or Yt- For reasons to become clear, it is convenient to write 


x = {x,X) and ]l={h,Y) 

where x E Y ■= x denotes the energy coordinates and X = (Xi,... ,Xm) 
the positions of the particles. Writing Xt = {x^,Xt), we observe that Xt consists of 
M independent continuous-time random walks on Fl, with “reflection” at &Dl as 
deflned earlier. Likewise, nEy:= x x gives the number of packets at 
the various sites or carried by particles, while Y = (Yi ,..., Ym) denotes the positions 
of the particles; and if Yt = (n^,!*), then Yt has the same description as Xt. Notice 
that F is really a function of {x,n). 

We now formulate a version of pathwise duality, counting sample paths of Xt (or 
Yt) in the following way: We say a “move” in Xt occurs when either (a) a particle 
jumps from z E T>l to w E Vl or (h) it jumps - all in one instant - from z E Vl 
to tc E Bl and back, and we regard different w E Bl as corresponding to distinct 
sample paths. Now flx a time interval (0,r), and let ti < ^2 < • • ■ < be the times 
in (0, r) at which Xt moves. To avoid discussing tf (i.e. just before or after tt), we let 
So = 0, Sn = T, flx arbitrarily s* E (ti, U+i) for i = 1,..., n —1, and agree to abbreviate 
this sample path as cr = (no, Ui,..., (T„) where Uj = W., with the understanding that 
in case (b), both the particle and the bath location involved are specifled. We also 
use the notation cr~^ = (a^,..., Uo) to denote the sample path corresponding to cr 
parameterized backwards in time, i.e., for cr~^, moves are made at times ti < ■ ■ ■ <tn 
where ti = r — tn+i-i] and at times • • •, the moves of cr“^ are the reverse of 
those in cr at times F, tn-i, • • • • 

Our one-step duality lemma reads as follows: 


Lemma 6. For any fixed x E X, n E y, and any sample path a = 
we have 


E(F(h,h.,.)|hQ = h, cr) = E(F(h.,.,h)|hQ = h, cr . 


((Jo,ai) on [0,r], 


Proof. The sample path a = (no, cri) describes exactly one move on the time interval 
(0,r). We consider separately the two cases corresponding to the two terms in the 
generator G of Xt (see Sect. 2.1). 


Case 1. Particle k jumps from site x E T>l to site w E T>l- The term corresponding 
to k and w in GiF can be written as I ■ II, where 


I = 


LttILtt Y TT \t(Z 

2<iD njD; fi,! D [ \L 
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and 


II = f / (L + p"'(i - pf-dp- iT't 

\J0 


From this we deduce that 


1 

1 

nfi. 

hfc! 

1 

1 

nz\ 

hfc! 

1 

1 


pi 


{^z + Vk) 

^n^+nfc 


nz+fik 


nj Hkl 
1 




^"■^(1 — p)"''dp 

Ctk 


{jlz d' llfc)! Yiz+n)s—m 


n^ + Uk + l 


^ m\{nz + fik — m)\ 


rtzj^k ^rn nz+nk-m 

Pk 


m! (n^ + nfc - m)! 
m=0 ^ ^ 


np.Ukl 


{nz + Uk + 1)! 


Thus I ■ II is the term corresponding to AiF in the generator of Yt with indices k 
and z. 


Case 2. Particle k jumps from z G T>l to w ^ Bl and back. The term corresponding 
to k in G 2 F and w E Bl can again be written as I ■ II, where / is as above and 


II = 


1 1 
nj. hk\ 


uo 


dp[p& + %))”• Tdor*/? (I) -C'!)?* 


From this we obtain 


rr _L _ 

= (^ / dplp K, + rjt)]’ 


'0 


riz + l 






m 

z 'Ik 


m\ (uz — ni)\ 

m=0 '' ^ ’ 



In the last equality, we used the simplihed version of the computation of // in Case 
1 (by letting hk = 0) and computed an elementary integral. We conclude that I ■ II 
is the term corresponding to A 2 F for the indices k and z. □ 


Next we extend Lemma El to sample paths involving arbitrary numbers of moves. 

Lemma 7. For any fixed x G fit, n G y, and any sample path cr = (cxo, cxi,..., am) on 
[0, r], we have 

^{F{n,xfi}\xQ = X, a) = ^{F{n^,x)\nQ = n, cr~^) . 

Proof. We prove by induction on m, the number of moves. The case m = 1 is Lemma 
1 Assume that we have proved the statement for < m — 1 moves. Letting si be as 
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defined above, we have the following: 

E(F(n,h^)|ho = x, (cro,cri,... ,(7^)) 

= j E{F{n,x^)\x^^ = F, (di, = h'l^o = (^-Q, cri))dh' 

= [ E(F(n^_^^,x')|no = n, (am, ■ P(x,^ = h'|ho = (^-Q, cri))dx' 

F(n',h') • P(n^_s, = ^IRq = n, (am, ■■■, cxi)) • P(x^^ = F\xq = x, (ag, cri))dh'dn' 

= J E(F(^,x^J\xq = x,(ao,ai)) ■ P(n^_^^ = ^\ng = n,(am,--,(Ti))dn' 

= j E(F(n^,x)\n^_^^ = h', ai, do) ■ = n'|no = n, ((T,„, ..., cri))dn' 

= E(F(n^,x)|ho = n, (am,(Tm-i, 



□ 

Remarks. The dnality statement above is a little more involved than that for Markov 
processes with disjoint phase spaces (see e.g. Proposition 1.2 in |JK14] i. Here, the 
phase spaces intersect in the set of particle conhgurations, and duality is proved one 
sample path in particle movements at a time; that is why we call it pathwise duality. 
Note that it is necessary to use reversed paths in the dual process to guarantee that 
Lemma E] can be applied from one step to the next. Also note that the expression 
’’pathwise duality” has been used in in Chapter 4 of |,TK14] in a different context: 
the ’’strong pathwise duality” and the ’’conditional pathwise duality” as dehned in 
jJK14] imply the usual duality while our dehnition is a weakening of that. Finally, 
we mention that we will have to further weaken our concept of duality for the proof 
of the case of systems with positive density of particles, see Section [3 

Duality has been used to prove LTE for a number of situations; see |CGGR13] and 
the references therein. Many of the ideas above have their origins in |KMP82j , though 
modihcations are needed as energy is not carried by particles in the KMP model. A 
similar set of ideas was also used in |RY07j . which considers the same model as ours 
in one space dimension and with only one particle. 

4.3. Consequences of duality. Let h* G A and G (V be hxed. Integrating over 
all sample paths <t of Xt on [0, t], Lemma [7] together with the fact that P[cr] = P[cr~^] 
imply that 

(10) J E(F(n^,x^)\xg = x^,(t) P(dcr) = J E(F(n^,xJ\ng = n^,cr-^) P(dcr"^) . 

Letting f —)■ oo, the left and right sides of the equation above tend to the two sides 
of the formula below: 
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Lemma 8. For any fixed G y, we have 

(11) A') = 


vGBl 



where pn^ is the asymptotic distribution of Yt as t ^ oo averaged over all Y with 
uniform distribution, assuming that Yq = {n^,Y). 


That the left side of flTOj) converges to the left side of ffTTj) as f —)■ cx) follows from 
the fact that the distribution of converges to (Proposition 2). The limit on 
the right side clearly exists, since all packets are eventually deposited in Bl, resulting 
in the simplihed form of F{ji,Yf) in the integrand. Notice also that the right side 
does not depend on consistent with the fact that the convergence to on the 
left is independent of initial condition. 

We now identify the relevant choices of In the proof of LTE for site energies, 
for example, we fix x G P, S' C and nonnegative integers (n*)yQS- If is chosen 
so that 

(12) = I if n ^ fiy = 0 Vv, and hj = 0 Vj , 


then the left side of flTT]) is equal to 


(13) 



^ (xL)-\-v 




a constant times the (n*)-moments of the distribution pj^^^ defined in Sect. 2.2. 

Thus the key to understanding p^^l is pn^. To get a handle on this distribution, we 
hnd that instead of working with Yt, which describes the evolution of the density (or 
distribution) of packets, it is productive to switch to an equivalent model that focuses 
directly on the movements of individual packets. Moreover, since only asymptotic 
distributions matter, we may work with a discrete time model, as long as the order 
of the steps are preserved. 


Discrete-time version of Yt focusing on movements of packets 

Consider a Markov chain Zk = fc = 0,1, 2,..., the variables of which are 

Zk = {Zi^k, • • •, ZM,k'-, Yi^ki ■ ■ ■ 1 YM,k) G {T>l I±) l±) {1, 2,..., M})^ x , 

with some fixed positive integer N (to be specified later) and the notation l±) for 
disjoint union. The hrst N coordinates of Zk describe the positions of the N (named) 
packets in the system, the position of packet i at step k being Zt^k, and the hnal M 
coordinates give the positions of the M particles (abusing notation slightly by using 
Y in both the continuous and discrete time models). The meaning of Zt^k ^ is 

obvious, and Zt^k = j means packet i is carried by particle j at time k. The transition 
probabilities of Zk are as follows: Given Zk, we choose with equal probability one of 
the M particles, say particle j, and choose with equal probability one of particle j’s 
neighboring sites, say w. If w G Vl, then we set Yj^k+i = it, and mix the packets 
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carried by particle j with those at site w by pooling them together and designating 
a random fraction of them to be carried by particle j and the rest to be left at 
site w. If ic G Sl, then all the packets carried by particle j are dropped off at 
site w, and Since all the packets are eventually dropped off in Bl, 

Zi^oo '■= limfc^oo Zi^k exists for 1 < i < almost surely. 

Let X E V, S G and L S> 1 be hxed. Associated with each 


n* = ■ ■ ■ ,n*M) 

is the Markov chain Zk whose initial condition Zq = (hg, Yq) is given by the following: 
hg is prescribed by n^, i.e. at time 0, there are n* packets at site {xL) + v and h* 
packets at site v E Bl, n* packets are carried by particle j, and Yq is uniformly 
distributed among all particle conhgurations in We claim - and leave it to the 
reader to check - that Yt and Zk differ only in the identity of individual packets and 
time changes that preserve the order of the moves, so they have the same asymptotic 
distribution, i.e. 


(14) 






where E is with respect to the evolution of the process Zk and 

M 

vGS v&Bl i=i 

For future reference. 


(15) E(-) = J E(^^,y)(-) f/(dF) 

where U is the uniform distribution over all particle configurations Y. 
We are now ready to prove Theorem 1. 


Proof of Theorem 1. Let be such that 


Uv 


1 if T = {xL) 
0 otherwise 


and fly 


0 for all n, = 0 for all j . 


Then the left side of CHi is equal to £^'^^(■^(3,^)), and the right side is given by 
J T (l;) dpn^, where pn^ is the asymptotic distribution of Zi^^o in the Markov chain 
above with N = 1 and Zi q = (xL). From the transition probabilities of Zk, it is 
clear that if we (i) disregard waiting times, i.e. times at which Zi^k does not change, 
and (ii) view the location of the packet when it is carried by particle j as Yj^k, then 
the trajectories of are those of a SSRW on V^. By Propositional as L —>■ cx) 
the distribution of z[^^ rescaled back to dV is the hitting probability of Brownian 
motion starting from x E V. Hence / T dp^^ —)■ u{x) where u is the solution of 
Laplace’s Equation with boundary condition T. □ 


Next we observe that Theorem 2 is reduced to the following proposition, the proof 
of which will occupy the next two sections. 
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Proposition 4. Let d >1, V C T on dV, x eV and S (ZU^ be prescribed. We 
fix also = {{nl)^^Si 0, 0,..., 0), and let be the Markov chain associated with 
n^. Then letting N = have 



Proof of Theorem 2 assuming Proposition We first prove the result assuming tight¬ 
ness of the sequence • Let /ioo be a weak limit point. Then putting 

together ffTTj) . flTT|) and Proposition 01 we see that the moments of /ioo are those of 
a product of exponential distributions with parameter fi{x) = u{x)~^. Hence /ioo is 
such a product; see Sect. 3.2. Since this is true for all limit points of /i^^ 5 , we conclude 
that /i^^g converges weakly to the measure claimed. 

It suffices to prove tightness one coordinate at a time, so we may assume S = {u}. 
Then the same reasoning (with n* = 2 in Proposition0]) implies that sup^, i^‘{xL)+v) — 

C for some C < oo. Chebyshev’s inequality then gives > n) < ^ for all 

L, proving tightness. □ 


We close this section with the following lemma. 

Lemma 9. Let be a system with N packets, Zfi be another system with iV -|- 1 
packets, and suppose both have the same number of particles. Assume further that 
Zf^Q = Zifl for i = 1,..., and = Yj^ for all j. Then with Zk coupled to 
the corresponding coordinates in Z^ in the natural way, we have Zf^, = Zi^^ for all 
i = 1,... ,N and k > 1. 


Proof. Without loss of generality, suppose that at step k, particle 1 jumps from site 2 ; 
to site w, and that the union of the packets carried by this particle or at site w prior 
to the mixing are labelled {1,..., n}. Then the probability that after the mixing, the 
set of packets carried by particle 1 is exactly {ji, C {1,2, ...,?7,} is given by 

11 lUn-iy. 

p{ji, ■■■iji) — —tttw “ 7—nrr 
n + l [n + 1)! 


If packet -|- 1 is not at site w, then clearly the situation is not disturbed. If it is 
there, we compute 


+PUl,-;jl,N +1) 


11^11 
n + 2 ("/■) + n + 2 («) 


l\{n-l)\ 
{n + 1)! 


Hence the dynamics of the first N packets are unaffected. They are also clearly 
unaffected if particle 1 drops off its packets at the bath. □ 


Remark. An implication of Lemma 0 is that when N > 1, the motion of each 
individual packet, when seen in the light of (i) and (ii) in the proof of Theorem 1, is a 
SSRW. Thus Proposition 0] is proved if these SSRW are independent, or close enough 
to being independent. This is what we will show. 
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5. LTE FOR SITE ENERGIES: d>2 

In Sect. 5.1, we introduce, mostly for convenience, a small modification of the 
process This modihed process is used a great deal in the pages to follow. Sect. 
5.2 contains the proof of LTE for site energies (Theorem 2) for d = 2. Due to the 
transience of SSRW, proofs for d > 2 are simpler and are given in Sect. 5.3, along 
with the proof of Theorem 4. 

5.1. A slightly modified process. We have seen in the proof of Theorem 1 that 
with a suitable modihcation of Zn, the movement of the packet becomes a SSRW. 
We now carry out the same type of modification systematically under more general 
conditions: 

The phase space of Zn is (VlUBl)^ , and its dynamics are derived from those for Zn 
in the following way: First, let = 7r(Zj^„) where n{Zi^n) = Zi^n if € T^l UI^l, 
and n{Zi^n) = ye,n if Zi^n = We then let to = 0, and for j = 1, 2,..., define 

tj = min e {1, 2,..., A} : Z' ^ Z' J . 

Finally, set Zj^„ = Z'^^. That is to say, hrst we confuse being at a site and being 
carried by a particle at that site, and then we collapse the times when there is no 
action according to this way of bookkeeping. 

Remark. We recognize that Z„ is not Markovian (and is not especially nice as a 
stochastic process). However, the order in which the N packets move about on 
VlU Bl is preserved as we go from Z'^ to Z„, even as time has been reparametrized. 
As a consequence, Z„ has the following important properties: 

1. For each L, the joint asymptotic distribution of (zj'^,..., Z^^) is identical 

2. Each packet individually performs a SSRW on Vl U Bl modulo waiting times 
(during which it stands still). 

3. The addition of new packets in the sense of Lemma 9 does not affect the order 
of movements of packets already under consideration. 

When two packets are at the same site, their next moves are not independent. We 
prove a uniform bound on how long they are likely to stick together: 

Lemma 10. Assume N = 2, and Zi^^o = -^ 2,^:0 ^ ^l- Let k he the smallest positive 
integer such that Zi^^^+k 7 ^ ■^ 2 ,fco+K- Then 

> k) < (2/3)^ A; = 1 , 2 ,.... 

Proof. The only way to hnd out what happens in the Z-process is to go back to the 
corresponding step in the Z-process. Below we enumerate all possible scenarios for 
Z‘i,tkQ)i = 1)2, that correspond to Zi^^o = ■^ 2 ,fco) and consider for each scenario the 
probability of the two packets staying together in the next one or two steps: 
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Scenario 1 . At time ko, exactly one of the packets is carried by a particle, or the two 
packets are carried by different particles. In both cases, a = 1, i.e., they will separate 
in the next step. 

Scenario 2. Both packets are carried by the same particle. Then P(fi; = 1) = 0 but 
P(a = 2) > i. Reason; This particle jumps, carrying both packets to the next site, 
where with probability | it drops one packet and carries the other, a scenario that is 
guaranteed to lead to a = 2 . 

Scenario 3. Neither packet is carried by a particle. Before the next move can occur, 
a particle has to enter the site, and with probability (|, |, |), picks up (i) neither, (ii) 
one, or (iii) both of the packets. If (i) occurs. Scenario 3 is repeated, (ii) and (iii) are 
followed by Scenarios 1 and 2 respectively. □ 

Notation: In this paper, we denote every universal constant by C, so that each 
occurrence of C may stand for a different number, even in the same line. 

5.2. Proof of Theorem 2: d = 2. We now focus on the planar case. Let Zn = 
be as in Proposition 4, and Zn the modihcation of Zn as dehned in Sect. 5.1. For 
/3,6 E (0,1) and i,j G {1,..., N}, we dehne 

Tij = min{n : dist(Zi,„, 

and Ti = min{n : dist(.Zj^„, (xL)) > 

In general, the dehnitions of Tjj and 71 depend on packets other than i and j (due to 
the way we collapse time when going from Zn to Zn), so let us first assume these are 
the only two packets present. 

Lemma 11. Consider V C and assume N = 2. Then for every 13 G (0,1) and 
S G (^^, 1 — /9) there is a constant C = C{f3, 6) such that 

C 

P(7i < Ti^ 2 ) < 7 Y 5 X for all sufficiently large L . 

Proof. Our plan is to write 

P (Ti < ri,2) < P (Ti < ri,2, ri,2 < + P (ri,2 > 

< P (Ti < + P (ri,2 > CL^+^-^) , 

and prove that each of the two terms above is < for large enough L. 

Consider first the second term. Decomposing the steps before ri ^2 according to 
whether .Zi „ = Z 2 ^n, we claim that 

(16) P(|{n : ^i,n ^ Z 2 ,n,n< n, 2 }\ > 

where |.| denotes the cardinality of a set. To see that, we let oq = 0, and dehne 
On = mm{k > Qn-i I Zi,fc f Z2,fc} and Un = Zi,a„ - Z2,a„ • 

Then the process Un is a planar SSRW as long as it is away from the origin. Whenever 
the SSRW would reach the origin, Un performs two steps of the SSRW thus avoids the 
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origin (more precisely, for any x and any e*, Cj G {(0,1), (0, —1), (1, 0), (—1, 0)}, 
P(f/n+i = x + ei\Un = x) = 1/4 if a; + ei 7^ (0,0) and F{Un+i = ej\Un = e^) = 1/16). 
Observe that by construction 

\{n I 7^ Z 2 ,„; n < ri^ajl = min{||f/„|| > L^} . 

n 

It is easy to show by the invariance principle that there is a p > 0 such that 

P(min{||t/„|| > L^} < > p 


holds for any starting position Uq with ||tlo|| < In fact, the left hand side converges 
to as L —)■ oo, but we will not need this. By induction, ||f/n|| reaches in time 
kL'^^ with probability at least 1 — (1 — p)^. Since 1 — d — 5 > 0, the choice k = 
gives dUD. 

Now it is well known that for two dimensional random walks, the number of returns 
to the origin up to time is O(logL) as L —)■ oo. Furthermore, formula (3.11) 

in jETbO] implies that the probability of the number of returns to the origin being 
bigger than log^ L is bounded by C/L^^^. Letting Ok be the duration that the 
two packets stick together at their fcth meeting, we note that these random variables 
are independent, and each is stochastically bounded by the geometric distribution in 
Lemma [TUI Thus we have 

P (|{n < ri,2 : = Z2,J| > 


< 


P 


meetings before ri ^2 > 


^100 

71 



c c 

£100 £100 ■ 




+ P 




m=l 


This completes the proof of P (ri ^2 > 

To hnish, we let Wn denote the SSRW corresponding to Zi^k with waiting times 
collapsed, and note that 


P (Ti < < 


< 


P| sup II VhFfcll > 

^<CLl+/3-^ ^ 


n=l 


k=l 


Since our assumption (1 — /5)/3 < 5 implies (3 + 5 > |(1 + /3 — 5), it follows by 
moderate deviation theorems such as Lemma [3] (by projecting onto one coordinate, 
for example) that each of the probabilities in the last line is bounded above by a term 
of the form < . Thus the sum is < 7 ^ for large L. 

The proof of the lemma is now complete. □ 


We now return to the case of N packets for arbitrary iV, and let % and be as 
dehned at the beginning of this subsection. 
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Lemma 12. Let D C Given any £ > 0, let (3 = 1 — 2e and 6 = e. Then for any 
i E {1,..., N}, the following holds for all sufficiently large L: 

(17) ^ for some n> % and some j ^ i) < 4{N — l)^ . 

Proof. First we claim that for each i,j E {1,..., A^} with i ^ j, we have 

( 18 ) < P,j) < > 

i.e., Lemma 11 in fact holds for any pair i,j in a process with N packets. To go from 
iV = 2 to general N, observe that while the dehnitions of 71 and Tij depend on the 
packets present (more time steps are collapsed when there are fewer packets), flTSD 
concerns only the relation between tij and 71, not the actual values of these random 
variables, and relations of this type are not affected by the presence of other packets 
as noted in the Remark in Sect. 5.1. 

Now for hxed j ^ i. Lemma 0] applied to Wn = ■^j,n — tells us that with 
probability > 1 — 3e, packets i and j do not meet after r, ^ for L large. (The result in 
Lemma 4 is not affected by waiting times.) Thus picking L large enough so the right 
side of flTS]) is < e, we have 

P(Zi_„ = Zj^ri e T^l for some n > %) 

< ¥{Ti < Tij) + P(Zj_n = Zj,n for some n > > Tij) 

<C £ T 3£ = As . 

Summing over all j ^ i gives flTT)) . □ 

For d = 2, we now prove Proposition 4, from which Theorem 2 follows as explained 
in Sect. 4.3. 


Proof of Proposition f (d = 2): We will prove 


(19) 




-» 0 


as L —)■ cx) 


inducting on N. The case of iV = 1 is Theorem 1. We now assume flT^ has been 
proved for a process with N — 1 packets, and note that when embedded in a process 
with N packets, the same asymptotic distribution holds for packets 1, 2,..., iV — 1. 

Let an arbitrarily small e > 0 be hxed. We consider L large enough for flTT)) to 
hold, and for such L, we dehne the stopping time S to be the smallest n > Tn such 
that Zisf^n = Zj^n £ for some j ^ N (setting S = oo otherwise), and dehne 
another process Zf^^ with the property that Zf^^ = Z^^n for n < S, and it is a 
SSRW independent of the movements of the other — 1 packets after time S. This 
ensures (i) Zff^ is independent of the joint distribution of Zj^oo,j = 1 ,..., A^ — 1 , 
and (ii) ^ ■^Ar,oo) < 4(A^ — 1 )^ for all large enough L. 
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Simplifying notation by dropping the superscript and abbreviating T as 

T{zj) and T as T{z]^), we have 

\E{Ul,T{z,))-u{xf\ < mufs,^T{z,))-{T{zM)-T{z%)))\ 

+ |E(nfj-/T(;.,)) ■ (E(T(;.;^)) -E(T(;.^)))| 

+ \E{Ul-^T{z,)) ■ (EiTiz^)) - u(a;))| 

+ |E(nfj-,'T(;^,-)) ■ - u{xf\ • 

Of the 4 lines on the right side, the hrst and the second = 0{e), due to the fact that 
T is bounded and F{T{zn) 7 ^ Tlz"^)) = 0(e); the third tends to 0 as L ^ cxd by 
Theorem 1 and Proposition 3; and the fourth tends to 0 as L ^ cxd by our inductive 
hypothesis. This completes the inductive step and the proposition. □ 

5.3. Related proofs. Details aside, the proof of the d = 2 case of Theorem 2 can be 
summarized as follows: We £x two distinct length scales, and L^, and consider 
Til fhe time it takes packet i to attain a net displacement of and Tjj, the time 

it takes packets i and j to separate by a distance of . Notice that „ is a 

SSRW except for the fact that the packets tend to stick together for a random time 
with hnite expectation when they meet. We then showed that 

(a) with high probability, % > Tij, due to the difference in length scale and 
also to the fact that O(logL), the number of encounters between packets 
i and j before Tij, is insignificant; 

(b) two packets that are apart are not likely to meet. 

It follows from (a) and (b) that after time T, the trajectories of packets i and j are 
effectively independent, and the desired result follows from the harmonic measure 
characteristization of hitting probability starting from x eT. 

The proofs below follow the same argument, with some simplihcations. 

Proof of Proposition 4 (d > 2). Transience of SSRW in d > 2 simplihes the estimates. 
Specihcally, let Wn be a SSRW in d > 2. Then the number of encounters in (a) 
can be estimated by the fact that given e > 0, there exists K such that 

encounters > K) < e , 

and (b) follows from the fact that given e > 0 , there exists Ci such that 

= 0 for some k\WQ > Ci) < e . 

Indeed the two length scales in the proof in Sect. 5.2 can, if one so chooses, be 
replaced by two suitably related constants Ci < € 2 - □ 

Proof of Theorem 4 for d>2. Theorem 4 differs from Theorem 2 in that the initial 
locations of the packets may be 0{L^) apart. As can be seen from the sketch of 
proof above, the following two places in the argument may be affected: (i) With 
regard to the number of encounters before Tij in (a), the probability of meeting at 



LTE FOR CERTAIN STOCHASTIC MODELS 


23 


least once cannot be increased if the packets are farther apart, and once they meet, 
the probability of meeting again is independent of their initial separation, (ii) For 
each L, when rescaled back to T> the packets do not start from x bnt from Xj with 
\xj — x| = The convergence of hitting probabilities starting from these 

slightly perturbed initial conditions is covered by Proposition [3l □ 


6. LTE FOR SITE ENERGIES: d=l 

In dimension 1, independent random walkers meet too often for the type of argu¬ 
ment in Section 5 to work. On the other hand, when two packets meet, it suggests 
the possibility of exchangeability, and we will make use of this in our proof. We may 
assume without loss of generality that V = [0,1], so that Bl = {—1, -h +1}, and write 
C = -l,7^ = L + l. 

6.1. The case of iV = 2. The argument in this subsection is borrowed from [RY07] : 
we need only the N = 2 case, which is considerably simpler. Starting from usual 
initial conditions, we let Ai = C or 71, and dehne 

Pl(AuA,) = P(zW =>1„Z<« =-4), 

Pi(T.*) = m[%=A,), p{;-) = lim etc. 

L^oo 

Lemma 13. It is sufficient to show that 

P{7Z,7Z) =x^ . 

Proof. We know from Id SSRW (or the gambler’s ruin problem) that 

(20) P{*,71) = X and P{C,*) = l — x. 

If P{71,71) = then P{C, 71) = P{*, TZ) — P{71,7Z) = x — x"^ = x{l — x), which is 
also equal to P{71, C) by symmetry. Thus as L —)■ cx), Zj^oo and Z 2 ,oo are independent 
Bernoulli random variables with weights (1 — x,x), giving 

□ 


Lemma 14. P{7Z,7Z) > x‘^ 

Proof. Let f{x,y) = xy. For each zi^^-process as dehned in Sect. 5.1, dehne 

^/(Zl,n5 '^2,n) ■= lE(/(Zi^„_|_i, Z2,n+l) | '^2,n) ~ f ^2,n) ■ 

Then A/ = 0 for all situations except when Zi_„ = Z 2 ,n, in which case we have 

E(A/ I two packets move together ) = ^((i — 1)^ + {i + 1)^) — = 1, 

where i is the location of the packets; otherwise A/ = 0. 
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Therefore, for any finite L, /(Zi_„, Z 2 ,n) is a bonnded snbmartingale. As n oo, 
f(Zi^n,Z 2 ,n) converges to /(Zi^oo, Z 2 ,oo) almost snrely. Fnrthermore, by the snb¬ 
martingale property, we have £'(/(Zi^oo, Z 2 ,oo)) > -E(/(Zi,o, Z 2 ,o))- Thus 

P(7^,7^) = lim L- 2 E(/(Zi,oe,Z 2 ,oo)) > lim L-^ ■ 

L^oo L^oo 

□ 

Lemma 15. P{TZ,TZ) < 

The proof of this lemma is simpler if we work with a process that differs from Z„ 
by a half-step. More precisely, in a packet first jumps before it picks up a random 
fraction of the packets at its destination. We let be with the order of jumping 
and mixing reversed, and apply the procedure at the beginning of Sect. 5.1 to to 
obtain a process we call Z„ (instead of Z„). Clearly, as L —)■ cxo, Z„, Z„ and all 
have the same asymptotic packet distributions. The notation P{■■,■) below refers to 
Zn- 

Proof of Lemma 15. We consider the function 

f{x,y) = xy- c\x - y\ 

where c is a constant to be determined. Let A/ be as before. Then it is easy to 
see that A/ = 0 when Zi^„ 7 ^ Z 2 ,n- If Zi^„ = Z 2 ,n, then A/ = 1 if the two packets 
move together in the next step, and A/ = —c if the next step involves exactly one of 
the packets. Now the two packets can move together only if they are both available 
to be picked up at step n -|- 1 , and even in that case, the probability that exactly 
one of them is picked up by the particle in question is = |. Thus we conclude that 
independently of what is going on in the Z^^-process, 

( 21 ) P(Zi^„+l 7 ^ Z 2 ,n+l|Zi^„ = Z 2 ,n) A ^ • 

Choosing c > 3 therefore will ensure that A/ < 0 at each step. Arguing as above, we 
then obtain P(7l,7l) < x^. □ 

Notice that the left side of can be zero if we use Z„ instead of Z„, for two 
packets carried by the same particle will necessarily move together if the particle hrst 
jumps before it mixes. This is our reason for using Z„. 

6.2. Proofs of Theorems 2 and 4 in the case d = 1. Following [KMP82] . we use 

the method of exchangeable random variables to extend the results above to the case 
of N packets for arbitrary N. 

Fix X G [0,1] and S' C Z. Let ^ 1 ,^ 2 , • • • G S (repeats allowed). For each N and 
L, we consider the process Zn = zlt^ with = (xL) + Vi,i = 1, 2,..., A, and for 
Aj = £ or 7?., we dehne 

■ ■ ■ 5 ^n) = P ^ = Ai, . . . , , 
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where the probability is dehned to be the average over all initial particle conhgurations 
Y. We begin with a lemma that sets the stage for exchangeability. 

Lemma 16. For every N E and every permutation a 0 / {1,, N}, 

lim (^ 1 ) • • • ) ^7v) PL,vi,...,Vff ) • • • ) ^(t(A) )) 0- 

L^OO ' V / V / / 

Proof. Since binary transformations generate the symmetric group, it suffices to con¬ 
sider a with the property that for some i ^ j, a{i) = j, (j{j) = i, and (j{£) = £ for all 
£ 7 ^ i,j. For hxed L, we consider the process Zn = zi!"\ and begin with the following 
observation; 

Abbreviating Pl,vi,...,vn p, we claim that 

p{Ai, ...,An\Y,F)= p{A^(i),A^(^ n)\Y, F) 
if F is one of the following two types of events: 

(i) F = {Zi^n = Zj^n} for some n; 

(ii) for two particles k ^ k' and for some n, 

F t^Z^Yi k^ Zj^n k , Y]^yi Yj^f and Z^^y^ ^ k^k 

To see that the asymptotic distributions are as claimed in case (i), we simply switch 

the roles of packets i and j from time n on. In case (ii), packets i and j are carried by 
two different particles, which are at the same site at time n. In this case, we switch 
not only the roles of packets i and j but the sets of randomness for particles k and 
k' from this time on. The condition that these particles do not carry other packets 
at time n ensures that the asymptotic distributions of other packets are not affected. 
We will refer to an event corresponding to (i) or (ii) above as a “viable switching” for 
packets i and j. 

Let Zn be the process obtained from Zn- As explained in Section 5, module waiting 
times, Zi^n—Zj^n is a SSRW when it is 7 ^ 0. When Zi^n = -^yn, the duration they spend 
together following each encounter is controlled by Lemma 10. As these durations are 
independent for different encounters and are bounded by random variables with hnite 
expectations, it follows from the discussion in Sect. 3.1 that the number of times 
packets i and j meet before reaching the bath is 0{L)-, we in fact need only that 

( 22 ) of encounters > C) ^ 1 as L —)■ cxd for every C. 

Let fi < f 2 < ... be the times packets i and j meet in in the Z^-process, i.e. 
Zi^r^ = Zj^r^ E T>l, and let Tq denote the corresponding times in the Z„-process. We 
let Fq denote the event that a viable switching occurs at the gth meeting, and will 
show that there exists b = b{N) > 0 such that 

(23) 1P(-^(?| ■ ■ ■, and Tq occurs ) > b . 

Once we prove this, the assertion in the lemma will follow: We make the relevant 
switch the hrst time a viable switching occurs, and the probability that this occurs 
before the packets reach the baths tends to 1 as L —)■ cxd. 
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To prove ([23]), suppose the gth meeting takes place at site v. Confusing with 
Zn momentarily, we assume for definiteness that packet i is the hrst to arrive at site 
V, where it remains through time Tg at which time packet j is brought to site v by 
particle k'. Under these assumptions, there are the following possibilities: 

Case 1. Zi^T-^ = V. Since = u) = |, P(-Fg) = | in this case. 

Case 2. = p, and = k'. In order for a viable switching to occur, packets 

i and j must be the only packets carried by particle k' at the end of step Tg. There 
being a maximum of N packets available to be picked up by particle k' at this time, 
we have 

(24) nZi,r, = Zj,,, = k' and Z,,,, ^ k'^ i,j) > , 


Case 3. Zi^r^-i = k for some particle k. Since k ^ k', for a viable switching to occur, 
packet i must be the only packet picked up by particle k the last time it moved; that 
probability is > jvp^qriy- As packet j must also be the only packet carried by particle 

/ \ 2 


k' at the end of step Tg, we have P(Uq) > 


1 


N(N+l) 


We have thus a lower bound b for Fg that depends only on N. 


□ 


Notice that the argument above applies equally well to the setting of Theorem 4: 
even with their initial positions 0{L^) apart, we still have fl2^ beacuse of the gambler 
ruin’s estimate and the argument above. 


Proof of Theorems 2 and 4 for d = 1. Let Vi,V 2 , - ■ ■ G S and Pl,vi,...,vn de as defined 
at the beginning of this subsection. For each N, Pl,vi,...,vn ^ convergent subse¬ 
quence as L —)■ cx) by compactness. It follows that there exists —)■ cx) so that as 

n —)■ CX2, 

Poc,v^,...,VM for every N . 

Equivalently, if x[^\X 2 ^\ ... are {0, l}-valued random variables defined in such a 
way that when {0,1} is identihed with {£, TZ}, the joint distribution of {x[^\ ..., X^'^) 
is equal to Pl,vi,...,vn fo^ ^^-ch N, then there is a sequence of random variables (Xi, X 2 ,...) 
to which {x[^ \ X 2 ^"\ ...) converges. 

By Lemmaiini the sequence (Xi, X 2 ,...) is exchangeable, and de Finetti’s Theorem 
applies (see Sect. 3.1). It follows from Sect. 6.1 that if m is the probability in Theorem 
5, then 

X = P(Xi = l)= [ 9dm , 

Jo 

= P(Xi = X 2 = 1) = f e^dm . 

Jo 

From these two integrals, together with Jensen’s Inequality, it follows that m = 6^, 
the delta function at x. That in turn implies, by the characterization of the measure 
m in de Finetti’s Theorem, that Xi,...,Xjv are independent. Since the analysis 
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above applies to any limit point of Pl,vi,...,vn as L — )■ cxo, we conclnde that Pl,vi,...,vn 
converges to the prodnct measnre as claimed. □ 


7. LTE FOR PARTICLE NUMBERS AND ENERGIES 


This section is dedicated to the proof of Theorem [3l 

We start with the following simple observation: The condition M(L)—)■ a vol(P) 
implies that for any n G "Dl, 


particles at site v = K} 
fM(L)\ f I f 1 
\ K ) ^ 


as L ^ oo. 


A similar compntation involving finitely many sites proves statement (1) of Theorem 


We will give a proof of statement (2) of Theorem |3] in the case S = {0}. This 
assnmption is not necessary, but it simplifies the notation considerably, and the proof 
of the general case is entirely analogous. 

To fix notation, let K G Z"*" be the number of particles at site (xL), and fix arbitrary 
nonnegative integers Uq, to be used as moments of the site and particle 

energies. For L with M{L) > K, let (7i,..., Tk) be an ordered list of distinct elements 
of {1, 2,..., M(L)}. We let Q(Ji, ...,Tk) denote the event that these are exactly the 
particles at site (xL), and define = h^(7i,..., Tk) by 


(25) 


= 


Un 


0 


if n = (xL) 
if V ^ (xL), 


n, 


= 0 Vn, 


Uj = 


n. 


if 7 = 7- 
otherwise . 


Lemma 17. The following holds for all large enough L: Let (Tl, ..., Tk) be fixed. We 
define as above, and let N = nQ + X)j=i ^\j] ■ Then 


(26) 




N 


Observe that the quantities on both sides of fl26|) are independent of (7i,..., Tr-), 
as any two lists of K particle names are clearly interchangeable, so the can be 
interpreted without ambiguity as convergence as L —)■ cx). Let ns denote the right 
side of (12^ by . [■^(a:)]^. 


Proof of Theorem 3(2) assuming Lemma [77| and S = {0}. Let 


= {# particles at site (xL) = K} . 


We project the conditional measure of pT on Ak, to the site and particle 

energy coordinates on (xL). The resulting probability, will be viewed as a 

measure on with coordinates (^,a;i, ...,uk), the .^-coordinate corresponding to 

site energy. 
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Let rig, be fixed. Using the notation above, we have 


^ u, 

11 ri* ' 

^0- =1 \y 


— y 

K\ ^ 


c'^o K [j] 

IIu^ ^Q(Tu-rrK)^{y\^y)\^K) 


{Ti,...,Tk)' 


Tlnl nj.A 

0 j=i b] 


y y -^(h*(7I, ...,Tff),i)lQ(ri,...,rif) 


(^r^^)-yiQK)[u{x)f--jj:^ = Kx)] 


{Ti,...,Tk) 

y T) 


N 


X)) 




The first two equalities are by definition. The convergence as L —)■ oo and the last 
equality are from Lemma [IT] and the comments following that lemma. 

As this holds for all rig, ...,nj'^j, we conclude, as in Sect. 4.3, that tends 
to a product of iL + 1 exponential distributions each with mean u{x). □ 


Proof of Lemma [13 Let (7i,..., Tk) be fixed. We say a particle configuration a E Q 
if QiTi, ..., Tk) holds, and introduce the function 

F'{n,x,a) = F(h,£)l{<TeQ} • 

The assertion in Lemma [T7| can then be rewritten as 

(27) j ~ e Q} ■ [n(a;)]'^. 


As in SectionjH we approximate the left side of fl27|l by x^, Xt)), and decom¬ 

pose into sample paths cr = (cxo, ..., am) of particle movements on [0, f\: 


(28) j E{F'{n^,x^,Xt)\^ = x^,(T)F{d(T) 

= y y HF{ru,Xt)\^ = x^,(T)¥{( 7 ) 

m cr=(ao,...,crm) 
o-m&Q 

(29) = y y E(F(hi,xJ|ng =n,,cr-^)P(o--^) . 

m (T=(ao,...,am) 

The first equality above is the definition of F' and the second follows from Lemma 
in Informally, for an X^trajectory that ends in a state with a set of particles at a 
certain site, its “dual trajectory”, which is obtained by reversing the paths of the 
particles, should start with the same set of particles at the same site. 
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Taking the limit t ^ oo in lines fl28l) and fl29|) we obtain 

(30) J F'(n„x,X)dli“-\i:,X) = ® 1^1 

To handle the right side of fl30|) . we alter the dehnition of in Sect. 4.3 slightly 
by restricting Y in Zq to Q, everything else unchanged. Let Uq denote this new 
probability distribution of?, and define E'(-) = J UQ{dY) (cf flT^ in Sect. 

4.3), so that the right side of fl30|) is equal to 

To obtain fl2711 . it remains to check that Proposition 4 holds with E replaced by E': 
For d > 2, this is not an issue, since properties of Y do not appear in the proof. For 
Id, one needs to check that the switching arguments in Sect. 6.2 are not affected by 
the restriction of Y to Q, and that is true as well. □ 




Appendix: Proof of Proposition [3] 

Our proof uses the continuous mapping theorem, which can be stated as follows 
(see e.g. Theorem 5.1 in |B68] ). As before, weak convergence is denoted by “=^”, and 
for a mapping / and measure p, /*/i is the measure given by /*/i(A) = p(/“^(A)). 

Theorem 6 (Continuous mapping theorem). Let X and X' he separable metric 
spaces, and let P{X) denote the set of Borel probability measures on X. We con¬ 
sider a Borel measurable mapping / : A —)■ A' with discontinuity set Df, and let 
pin, L ^ such that (i) pin ^ pi as n ^ oo and (ii) pi{Df) = 0. Then 

f*Ln as n ^ oc. 

To fix some notation, for T G we let 

A = C([0,T],M‘^) 

be the set of continuous maps from [0, T] to endowed with the sup norm, making 
it a separable metric space. For a G we let 5“ be the standard Brownian motion 
starting from a G up to time T, and with a slight abuse of notation, we use B°' to 
denote also the corresponding measure on A. Let hF“ be the rescaled SSRW up to 
time [nTj starting from a G i.e. if Sk is a d-dimensional SSRW with Sq = (ay/n), 
then Wf{k/n) = ^ for /c = 0,1,..., [riTj, and fF“(t) is obtained by interpolating 
linearly between t = k/n and t = {k + l)/n. As with we use to denote also 
the corresponding measure on A. By the invariance principle, 

(31) Wn B^ as n ^ oo 

on any interval [0, T]. The convergence in fIST]) differs from that asserted in Proposition 
13] in that the latter is for paths that terminate not at a hxed time but upon reaching 
dV. 
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Proof of Proposition 0 Let x eP and £ > 0 be given. We write 

pi = {y e EV,\y-z\< 1}, 

and let % = %{e) be such that a Brownian motion starting from x reaches dP^ before 
time T with probability at least 1 — £/(2||T||oo)- For u E X = C{\f), T], M'^), we dehne 
t{uj) = min{t G [0, T] : uj{t) E dP} if such a t exists, = oo if it does not. Then we 
dehne / : X —)■ M by 

fhd = / T{uj{t)) ifr<T 
I maXygax,T(i/) if r = cx) . 

Lemma 18. B^{Df) = 0 


We hrst hnish the proof assuming the result of this lemma. Via a rigid translation, 
we may assume a: = 0 (so that xL E if" for all V). Then with Sq = xL, Sn/L in 
the proposition is W^ 2 {n/Lf) in the notation above. Observe that we are now in the 
setting of the continuous mapping theorem: W ^2 ^ as L ^ oo is condition (i) 
in Theorem El and the assertion in Lemma [TS] is condition (ii). Thus the theorem 
applies, and its conclusion together with our choice of T gives exactly ([S]) in the case 
S'o = xL. 

To prove the full statement of Proposition [3l we observe that if the result was false, 
there would be a sequence Xk E with Xk ^ x and a sequence Lk ^ oo such that 
if is the SSRW on with = XkLk, then 



where r is the smallest n such that E Bl^.- Such a scenario cannot occur: Since 
= (xk—x) + Wf 2 , it follows from fl^ that W ^2 B^ on [0, T], and the argument 

in the last paragraph with W ^2 in the place of W ^2 gives the opposite inequality. □ 


To complete the proof, it remains to show that the discontinuity set of / has zero 
Wiener measure. 


Proof of Lemma {TR, First, we identify the discontinuity set Dj. If t{uj) = oo, then 
the trajectory of uj up to time T is bounded away from dP, hence / is continuous at uj. 
If t{u) < T, then liminf^^/^i^ u(a;') > t{uj) for the same reason, but the corresponding 
limsup can be strictly greater than ripj) if the trajectory of uj does not cross to the 
other side of dP immediately following t{u). More precisely, we have deduced that 
Df = {t = %}\J E where 

E = {uj ■. t{uj) < T and 3rj = rj{uj) > 0 s.t. uj{{t, r + //)) C P} 
where P is the closure of D. 

Clearly, {r = T} has measure 0, so it suffices to show B^[E) = 0. 

Since harmonic measure is absolutely continuous, the set of uj for which uj{t) lies 
at a point at which dP is not differentiable has measure zero. Let uj be outside 
of this measure zero set, and fix an orthonormal basis {ei,..., e^} of such that ei 
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is the outward normal to dV at ca(r). Then there exist K > Q and a neighborhood 
U of u{t) in such that 

(32) yef) implies yi < K\\{y2, ...,ydW ■ 

Recall that by the strong Markov property, starting from the stopping time r 
is a Brownian motion. In particular, by projecting this Brownian motion, which we 
call B{t), to the line parallel to ei and to the hyperplane spanned by 62 ,..., 6 ^, we 
obtain two independent Brownian motions, Bi{t) and Bd-i{t)- Since Bi{t)/\/t and 
Bd-i{t)/y/t have standard normal distributions, 

P(|Ri(t)| < 0 and P(||Rd_i(t)|| > 0 

as t 0. Choosing j, 0 so that 

5^p(lA((„)l < . 5^P(ll4-i((„)ll > fi-'VVF) < 00 , 

n n 

it follows from the Borel-Cantelli lemma that \Bi{tn)\ > and ||.Brf_i(tn)|| < 
tn^/^/K hold for all but hnitely many n. For for which these inequalities hold, we 
are guaranteed that 13{tn) ii Bi{tn) > 0 and 13(tn) G U. 

Now it is a well known fact that on the time interval [0,p] for every p > 0, a 
ID Brownian motion starting from 0 makes inhnitely many excursions from 0, and 
each excursion is positive with probability 1/2 independently of other excursions. 
Applying this fact to Bi{t), and assuming (as we may) that each tn lies in a different 
excursion, it follows that with probability 1, 13i{tn) > 0 for inhnitely many n. Since 
P(.B(t) EU^t E [ 0 , 77 ]) ^ 1 as p —)■ 0, we have proved that B^{E) = 0. 

□ 
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